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Abstract 

We construct new classes of vortex-like solutions of the CP N model 
in (3+1) dimensions and discuss some of their properties. These solutions 
are obtained by generalizing to (3+1) dimensions the techniques well es- 
tablished for the two dimensional CP N models. We show that as the total 
energy of these solutions is infinite, they describe evolving vortices and anti- 
vortices with the energy density of some configurations varying in time. We 
also make some further observations about the dynamics of these vortices. 



1 Introduction 



In this paper we present new classes of vortex-like solutions of the CP model 
01 [2J in (3+1) dimensions. Our results generalize those obtained in our previous 
paper [3 J where we presented a quite large class of exact solutions of CP N models 
in (3 + 1) dimensions. These solutions were described by arbitrary functions of 
two variables, namely of the combinations x 1 + i x 2 and x 3 + x°, where x^, /x = 
0, 1, 2, 3 are the Cartesian coordinates of four dimensional Minkowski space-time. 
Then we considered field configurations, which for fixed values of x 3 + x° were 
holomorphic solutions of the CP N model in (2+0) dimensions. The dependence 
on x 3 +x° was assumed to be in terms of phase factors (e tk<yX +x >). These solutions 
then described straight vortices with waves traveling along them with the speed of 
light. Solutions of that type were also constructed for an extended version of the 
Skyrme-Faddeev model Our previous paper [3| contained other solutions 

for which the vortices and the waves were in more complicated interactions with 
each other. 

In this paper we generalize the procedure of and generate many more 
vortex-like solutions and also discuss solutions which correspond to configura- 
tions of parallel vortices and anti- vortices. Such structures interact with each other 
and our solutions describe this interaction and the resultant dynamics. A novelty 
of the paper is that we generalize to (3 + 1) dimensions a method for constructing 
solutions which was originally proposed [2J in the context of the two dimensional 
CP N model. Given an holomorphic solution, i.e. a configuration depending only 
on x 1 + i x 2 and x 3 + x°, we are able to generate, using a projection operator, 
solutions depending on x 1 + i x 2 , x 3 + x° and also x 1 — i x 2 . 

As our solutions describe vortices their total energy is infinite so to compare 
various configurations of vortices it is convenient to talk of energy density or 
energy per unit length. Then, as we discuss in this paper interesting phenomena 
that can take place - the energy per unit length can stay constant, be periodic in 
time or even grow with time. At first sight this may seem surprising but, in fact, 
this is not in contradiction of any principles, as the total energy remains infinite 
and so is "constant" (i.e. does not change). This observation complements the 
observation of our previous paper fl3j in which we pointed out that although the 
energy per unit length of various parallel vortex configurations can depend on the 
distance between them the vortices would still remain at rest. 

The paper is organized as follows. In the next section, for completeness, we 
introduce our notation and recall some basic properties of the CP N models and 
of their classical solutions in (2+0) dimensions. 
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The next section presents our solutions and the following one discusses some 
properties of these solutions. We finish the paper with a short section presenting 
our conclusions and further remarks. 



2 General remarks about the CP N model 

The CP N model in (3 + 1) dimensional Minkowski space-time is defined in terms 
of its Lagrangian density 

£ = M 2 (D tl Z) t D' x Z, Z ] -Z = l, (2.1) 

where M 2 is a constant with the dimension of mass, Z = (Z%, . . . , Zn+i) 

and it satisfies the constraint Z^ ■ Z = 1.. The covariant derivative acts on any 

iV component vector \I/ and so also on Z, according to 

= - (z+ • d^z)m. 

The index fi runs here over the set /i = {0, 1, 2, 3} and the Minkowski metric is 



(+,-,-,-). The Lagrangian (2.1 ) is invariant under the global transformation 

Z —7- U Z, with U being a (iV+ 1) x (JV+1) unitary matrix. One of the advantages 

of the Z parametrization is that it makes this U(N + 1) symmetry explicit (DEI. 

It is also convenient to use the 'un-normalized' vectors Z with components Zj. 

Then 

Z = -r^=, (2.2) 

'z\-z 



where the dot product involves the summation over all (N + 1) components of Z. 

Sometimes, exploiting the full projective space symmetry of the model, we set 
u = - and so use the parametrization 

Z = ^^'■■■'^ (2.3) 



1 + \ui\ 2 + . . . + \u N \ 

The M-field parametrization does not make the U(N + 1) symmetry explicit but it 
has the advantage that it brings out the real degrees of the freedom of the model. 



In terms of Uj's the Lagrangian density (2.1 ) takes the form 

AM 2 



{l+u ] -u)d^u ] -d^u~(d^u ] -u){u ] -d^u) . (2.4) 
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The classical solutions of the model are given by the N Euler-Lagrange equations 
which take the form: 

(1 + u ] ■ u) d^d^Uk - 2(u ] ■ d»u) d^u k = 0. (2.5) 

The simplest CP 1 case is given by one function u: Z = j 1 '"- 1 . 

V 1+ M 2 



3 Some solutions 

In this paper we shall use the notation of i.e. we define 

z = x 1 + ieix 2 , z = x 1 —ieix 2 , y± = x 3 ±€2X° (3.6) 

with e a = ±1, a = 1, 2. 

It is easy to check that any set of functions Z k and so uu that depend on coor- 
dinates x^ in a special way, namely 

u k = u k (z,y+) (3.7) 

is a solution of the system of equations ( |2.5 ). The Minkowski metric in the coor- 



dinates ( |3.6[ ) becomes ds 2 = —dz dz — dy + dy-. It then follows that (13/71) satisfies 
simultaneously d^d^Ui = and d^Uid^Uj = for all i, j — 1, . . . , N. Hence this 
class of solutions is quite large. 

However, these are not the only solutions we can construct very easily. In 
fact, we can exploit the construction of the solutions of the CP N model in 
(2+0) dimensions (for N > 1) to obtain further solutions. To do this we recall the 
construction in (2+0) dimensions: 

First we define a Gramm- Schmidt orthogonalising operator P z by its action on 
any vector / £ C N+1 , namely 

PJ = dj - f -jj^-- (3-8) 
Then, if we take f — f(z) and consider Z — f(z) the corresponding u solves 



the equations (2.5). Note that as f(z) does not depend on y± we have a solution 
of the CP N model in (2+0) and in (3+1) dimensions. However, as is well known, 
(see e.g. [|2l and the references therein) 

Z = PJ(z) (3.9) 



3 



defines further it's which also solve (|2.5[) in (2+0) dimensions. But, as the expres- 



sion for u does not depend on y± these functions also solve the equations (2.5 ) in 



(3+1) dimensions. This procedure can then be repeated, namely we can take 

Z = P k z f(z), (3.10) 

where P k J = P^P^f). 

To have more general solutions we observe that, like in [3], we can make the 
coefficients of z in the original f(z) to be functions of one of y±, say, y + . As y + is 
real the operation of applying P z operator does not introduce the other y±, i.e. ?/ , 
and so the corresponding Z and so u give us further solutions of the equations ( ]2.5[ ) 
in (3+1) dimensions. This way for iV > 1 we can have holomorphic solutions and 
also 'mixed' solutions. 

They are given, respectively, by 

u k {z,y + ) = - r (3.11) 

f N+1 (z,y+) 

and 

pi 

u k (z,z,y + ) = z k . (3.12) 

Note that hke in the (2+0) case the last (as we take larger /) nonvanishing 
solution would be antiholomorphic. Then the corresponding will be functions 
of only z and y + . 



3.1 Some properties of our solutions 

Let us first discuss briefly some quantities which we will use in the discussion of 
various properties of our solutions. 



3.1.1 The energy of the solutions 

The Hamiltonian density of the CP N model, when written in coordinates (z, z, 
y + , yJ), takes the form 

n = H w +H {2) , (3.13) 
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where 



8M 2 



(1 + wt - u ) 2 

8M 2 
(1 + ut-w) 2 



4« f • A 2 • d z u + d z u ] ■ A 2 ■ d z u] (3.14) 



<W • A 2 • d+u + 0_«t • A 2 • <9_ 



(3.15) 



and A 2 , = (1 + u ] ■ u)8. 



um 



r 



For solutions depending on y + i.e. described by u k (z, z, y + ) the part of the 
Hamiltonian density ( 3.15| ) that contains 9 drops out. For the holomorphic solu- 
tions the second part of (3.14[) also drops out. For the 'mixed' solutions described 



by (3.12) both parts of (3.14) are nonzero. 

Note that as our solutions depend on variables x° and x 3 only through the 
combination y + it is useful to define the concept of energy per unit length which 
involves the integration over x 1 and x 2 (i.e. over the plane perpendicular to the x 3 
axis). This gives us 



£= / dx 1 dx 2 'H = 87iM 2 



J(l) +J (2) 



where 



x (a) = — - / dx 1 dx 2 n (a) , 

8?rM 2 Jri 



a = 1,2. 



3.1.2 The topological charge 

As we are working with vortex configurations it is important to introduce the two- 
dimensional topological charge defined by the integral 



Q 



top 



dx 1 dx 2 p top 



(3.16) 



whose density is given by 



Ptop = -£ij(DiZy ■ (DjZ) 



1 <W • A 2 • di 



7T 



Si 
7T 



1+Ut. u )2 



1 d s u^ ■ A 2 • d z u — d z v) ■ A 2 • d z u 
7T (1 + id ■ u) 2 



(3.17) 



The indices i and j here only take two values {1,2}. It is easy to see that for the 
holomorphic solution Q top —X^ 1 '. 
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4 Vortex solutions of the CP N model and some of 
their properties 

In we studied some general classes of solutions of the CP 1 model. Here, first 
of all, we concentrate our attention on two classes of holomorphic solutions of 
the CP 1 model and then look in some detail at the CP 2 model concentrating our 



attention this time on 'mixed' solutions (3.12). 



4.1 CP 1 solutions 

In the CP 1 model we have two functions /i and f 2 and in our discussion we can 
take their ratio w = 

h 

Let us first consider the case when all the dependence on y + is in the form of 
phase factors e %klV+ where ki are constant. Many interesting features are observed 
for the configurations given by 

h(z,y + ) = z 2 + a lZ e ik ^\ f 2 (z,y + )=a 2 z + a 3 e lk2y +, (4.18) 

where we have assumed, for simplicity, that all three parameters a%, a 2 and a 3 are 
real. The generalization to their complex values does not bring anything new to 
the problem. 

The holomorphic solution u is then of the form 

z _|_ ai e ikiy+ 

u{z,y+) = z ^— • (4.19) 

a 2 z + a 3 e ll * 2V + 

The zeros of denominator do not lead to the singularities in the energy density 
as both integrals and are invariant with respect to the inversion u—t-. 



Next we look in detail at various special cases of this solution (4.19). 



4.1.1 The tube solution 

First we consider the case of a\ = a 2 = 0. In this case the field configuration 
becomes 

z 2 

u = —e- lk2y+ . (4.20) 

It is easy to convince oneself that this field configuration describes a vortex 
with waves traveling along it with the speed of light. The profile of the energy 
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density is independent of y + . It has a maximum at a ring of radius r which 
satisfies r\ < r < r 2 , where T\ = ^J~^ is me radius of the circle at which the 

Hamiltonian density has a maximum, and r 2 = \J\a^\ corresponds to the 

radius of the circle at which Ti^ has a maximum. The radius tq depends on as 
and k 2 . For fc 2 — >■ it tends to n and for k 2 — > ±00 it tends to r 2 . The integral 
describes the topological charge of the vortex which for the solution considered 
here is 

J(D = 1 / dx l dx 2 





Z 


2 


(«3 + 


z 


4)2 



71 JR? 

The contribution to the energy per unit length that comes from the traveling 
waves can be also calculated explicitly. We find 

iS ^ — f dx dx k>o ~, o i TTTTT ~~ Ajo I Q>^t I . 

71 Jr 2 (a| + \z\ 4 ) 2 4 21 1 

A modification of a solution of this type had been already studied in [3J. An 
example of such a solution is shown in Fig [TJ where we plot the components of 
the i so vector 

n = — (u + u*, -i(u - u*), \u\ 2 - l) (4.21) 

1 + \u\ 2 V ' 

which depend on y + . As y + changes the images in Fig{T]rotate. In Fig. [2} we plot 
the two contributions, topological and wave, of the energy density on the solution 



(4.20). 



4.1.2 The spiral solution 



A less trivial but still a very simple solution is obtained from ( 4.19| ) by putting 

a 3 = 0, and so u is given by 

u = —(z + a 1 e ikiy+ ). (4.22) 
a 2 

In this case the integrals ZW and can be calculated explicitly. They take 
the values 

Z« = 1 / dx l dx 2 \- — ]-— = 1, (4.23) 



7r Jr? a\ (1 + \u\ 2 ^ 2 



*2 

,21,2 



J(2) = I / dx'dx 2 ^ - aik ' = a\k\. (4.24) 
Ti Jr 2 a 2 (1 + \u\ 2 ) 2 
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In Fig. [3]we plot the components of the isovector (4.21 ) for the solution (4.22 ) 
In order to analyze the energy density let us introduce the parameterization z = 



re ltp . Then 



\u\ 



+ a\ — 2a x r cos(y? — kiy + — ir) 



We note that the energy per unit length (H integrated over the x 1 x 2 plane) does 
not depend on a 2 or y + , whereas the energy density H does. The maximum of 
the energy density (\u\ 2 = 0) is located at r = |ai| and cp = k±y + + it. The 
curve (a x cos (k 1 y + + 7r), a\ sin (kiy + + 7t),y+) that joins the points at which the 
energy density has a local maximum is a spiral. On this spiral not only H has a 
maximum but so do also both its contributions and T-L^ 2 \ As y + = x 3 + x°, 
we note that the spiral rotates around the x 3 axis with the speed of light. The 
only effect of the dependence on y + is the rotation of the energy density. Thus the 
energy per unit length calculated for e.g. y + = is also valid for other values of 
the variable y + . 



4.1.3 The general case (4.19 1 



For general values of a±, a 2 and a 3 the expressions for the contributions to the 
energy become rather complicated. We can write them as 



1 



.4 



dit)C didC y^-yri " 

7T JR 2 C 2 



JW = - [ dx l dx 2 ^-. 

IX JR 2 C 2 



(4.25) 



where the expressions for A, B and C take the form (written in cylindrical coor- 
dinates (r, (p, y + ) with z = re llfi ) 

= a\a\ + 4agr 2 + a\r A + 2aia 2 a 3 r 2 cos [2<p — (hi + k 2 )y+] 

+ Aaialr cos (ip — kiy + ) + Aa2a 3 r 3 cos (ip — k 2 y+) (4.26) 

= r 2 a\al{ki - k 2 f + r^{a\a\k\ + a\k\) 

— 2aia\{ki — k 2 )k 2 r 3 cos [<p — kiy + ] 

+ 2aja 2 a 3 (k 1 - k 2 )kir 3 cos [ip - k 2 y + ] 



A 



B 



C 



2a x a 2 a z k x k 2 r^ cos \{k\ - k 2 )y + ] 
r 2 + 2a\r cos [ip — kiy + ] + a\ 



+ 4 



r 2 + 2 — r cos Up — k 2 y + ] + 
a 2 



(4.27) 



(4.28) 



To fully analyse these expressions requires numerical work. In Figs. [5] and [6] 



we present the plots of the isovector (4.21 ) as well as of the energy densities for a 
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particular example of the above solution. However, even for a general configura- 
tion, it is possible to make a few analytical observations: 

• Rotations: 

Note that the energy per unit length depends on y + through periodic func- 
tions, involving four frequencies, namely k\, k 2 and k\ ± k 2 . However, one 
can isolate four situations where only one frequency is relevant and the time 
evolution reduces to a rotation around the ,x 3 -axis. In such cases, A, B and 
C depend on if and y+ only through the combination if — coy + , and the four 
possibilities when this happens are: 



1. 


ki 


= k 2 = k and oj 


2. 


a i 


= and oj = k 2 


3. 


a 2 


= and u = ki 


4. 


a 3 


= and u = ki 



k 



Note that the spiral solution (4.22) belongs to the last case and the tube so- 



lution ( 4.20| ) corresponds to the case when none of the frequencies matters. 
Singularity 



S3-. Indeed, 
2nn, with 



The solution (4.19) exhibits an interesting property when a± 
in this case it reduces to u = z/a 2 whenever [k 2 — ki)y+ 
n integer. The case k% = k 2 is not interesting since it leads to a solution 
independent of y + . However, for k\ ^ k 2 the solutions change their prop- 



erties, including the two dimensional topological charge (3.16), whenever 
y + = £ n = k ^™ 1 ■ For those special values of y + the quantities (4.26 )-(4.27 ) 
become 



A — a 2 \r—r n \ , 



B = aKh-hfr^f-fnl 2 , C = (r 2 +a 2 2 )\f-f r , 



where f and r n are two-component vectors: f — > (x, y), and f n — >■ (x n , y n ), 
with 



03 

a 2 



COS (k^ n 



7T 



«3 
Cl 2 



sm(k^ n + tt). 



(4.29) 



The expression |r 



I 2 then becomes 



r — r r , 



(x - x n ) +{y -y n 

«3 



r — 2 — r cos ((f — ki£ n — n) 



«2 



as\ 2 
a 2 J 



(4.30) 
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The cancelation changes the degree of polynomials of variable z which 
causes the topological charge to jump from Q top = 2 down to Q top = 1. 
The new topological charge is then given by the integral 

Qtop = Z« = - [ dx l dx 2 . a " = 1. (4.31) 

Of course, such behaviour is well known from the study of topological soli- 
tons [6J. The space of parameters of the field configuration is not complete 
(has 'holes') and the integrand of the charge density has corresponding delta 



functions, which are not seen in (4.31 ). The interesting property here is that 



this process of the vortex shrinking to the delta function and then expanding 
again is a function of time; i.e. is part of the dynamics of the system and is 
described by our solution. 

The second and related important fact comes from the study of the integral 
Z( 2 \ One can check that when the vortex shrinks to the delta function (i.e. 
the cancellation takes place) the integral 

tf» = i / tftf ,^. (4.32) 
7r Jr 2 [r 2 + a,2) \t — r n \ z 

diverges. This divergence comes from the singularity at the point r = 
r n which is responsible for the energy of the solution becoming infinite. 
Clearly, from a physical point of view such field configurations should be 
excluded. 

Anti-holomorphic solutions 



We can now also apply the transformation (3.8 ) to (4.19 ) and this would give 
us an anti-holomorphic solution. Its properties are not very different from 
what we had for the holomorphic one (except that the choice and meaning 
of parameters is different) so we do not discuss it here. 



4.1.4 Further Comments 

In our discussion so far we have assumed that all y + dependence of the 2-dimensional 
u(z) is of the form of phase factors exp(iky + Ys. There is, of course, no need to 
be so restrictive. We could make the parameters of the 2-dimensional u(z) depend 
on y + in a more general way. Thus we could consider, for instance, also 

u(z,y+) = A — r, (4.33) 

z - a(y+) 
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where a(y + ) is an arbitrary function. 

Then, taking e.g. a(y+) = ay + would result in a vortex located at x 2 = 0, 
x 1 = ax 3 moving in the x 3 direction with the velocity of light. Taking a more 
complicated function, e.g. a(y + ) = ay\ would result in a curved vortex x 1 = 
a(x 3 ) 2 etc. One can also combine this dependence, for systems of more vortices, 
with the other dependences discussed above. This complicates the discussion but 
does not change its main features, hence in the remainder of this paper we return 
to the discussion of the dependence on y + through the phase factors. 

One could naively think that infiniteness of the total energy of our solution is 
related to some "improper" choice of the dependence on y + . This is not true since 
the origin of the divergence comes from the topological nature of W^. The fact 
that 'HS 1 ' is a total derivative prevents the dependence of Ti.^ on any parameters 
(including any depending on y + ). One can note that for some special cases like 
u = z 2 exp (—ay+) the contribution to the total energy coming from is finite 
but the total energy remains infinite since Ti.^ contribution is always present. 

4.2 The CP 2 model 

Next we consider solutions of the CP 2 model. First we look at the holomorphic 
ones. 

4.2.1 The holomorphic solutions 



The simplest CP 2 model solution can be obtained by adding to the system (4.18) 
a constant third function, i.e. define 

f 1 (z,y+) = z 2 + a lZ e ik ^ 
f 2 (z,y+) = a 2 z + a 3 e ik2V+ 

f 3 (z,y + ) = a 4 . (4.34) 

Then we can define holomorphic configurations as Ui = -4, i = 1, 2, i.e. 

z 2 + a lZ e ik ^+ a 2 z + a z e ik ™+ 

u 1 {z,y + ) = , U2 (z,y + ) = . (4.35) 

Alternatively, we can interchange f 2 -H- ^3 and consider the holomorphic config- 
urations 

~ ( \ z 2 + a x ze ik ^+ _ a A 

Ui(z,y+) = ■ tj—-, u 2 {z,y + ) = — -. (4.36) 

a 2 z + 03 e lK2y + a 2 z + 0,3 e ltt2y + 
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Note from ( |2.3[ ) that such an interchange corresponds to a phase transformation 
in Z, so both configurations describe the same solution of the CP 2 model. Note 



also (easier from (4.36)) that when a 4 — > this CP solution reduces to the 



holomorphic CP 1 solution discussed before. In fact, u 2 vanishes, and u x becomes 
the CP 1 u-field. 

The integrals for this CP 2 holomorphic solution (using definition 



( |4.35[ ) or ( |4.36| )) now take the form 

A 



j(D 



71 JB? 



C 2 



J(2) 



dx 1 dx 2 



71 JR-2 



B 

C 2 



(4.37) 



where A, B, C differ from A, B, C given by ( |4726| ), ( 1427] ) and ((4728J) by terms 
proportional to a 2 , i.e. 

A = A + a 2 A [al + al + 4r 2 + 4a 1 r cos (ip-k 1 y + )} (4.38) 
B = B + a\[a\k 2 r 2 + ajk 2 ] (4.39) 
C = C + a^. (4.40) 

The Hamiltonian density T-i^, which is proportional to g|, is now regular at r = 

(where previously we had a singularity) as now it 



f n and y + = £ n for ai 
takes the value 



03 

(12 



B_ 

C 2 



f=f n ,y+=£,„ 



a 
a 2 



4^1 + ^2 



Hence we note that going to the CP 2 manifold (by taking a 4 7^ 0) has 'filled in 
the hole' in the space of parameters (i.e. as the system evolves none of its vortices 
shrinks to the delta function). 

Note also that the energy density is independent of y + in four cases: k\ — k%, 
a\ = 0, a 2 = and a 3 = 0. 



4.2.2 The mixed solution 



Next we look at the 'new' mixed solutions. First we use (3.8) to calculate P z f 



We find that for the system (4.34) they take the form 



Pzh 



Pzh 



ale ikiy+ 



2ze~ ikiy+ + 



a 1 



+ e 



a 3 + a 2 ze ik2V+ 



ai a 3 + 2a 3 ze- iklV+ + a 2 z 2 e~ i{kl+k2)v+ 



2 

a 2 a 4 



12 



e ik2V+ z 



a x + ze ikiy+ 



a\a 3 + 2a 3 ze 



-ikiy+ 



+ (I2Z e 



2 -i(k 1 +k 2 )y+ 



Pzh 



-a 4 e 



-ik 2 y + 



2 -^ifc 2 y+ _|_ ^ e ik 2 y + (^ e ik iy+ 



a 2 a 3 + ape + ze 



+ ai){2ze 



-ik\y+ 



When written in terms of U; t this mixed solution is given by 



P z fi 



u 2 {z,z,y + ) 



Pzh 

Pzh 



(4.41) 



Note that in the limit a 4 — > the mixed solution ( |4.41[ ) becomes the anti- 
holomorphic solution of the CP 1 model mentioned before. However for a 4 7^ 
the solution is different. This time the expressions for the energy density are quite 
complicated - so we do not present them here. However, we note that to guarantee 
the convergence of the integral 1^ we have to require that a 2 7^ 0. 

To demonstrate that the energy per unit length does not depend on y + can be 
checked without much effort. First, we observe that the overall factors e %k]V + in 
P z fk do not matter as they cancel in the expressions for \uj\ 2 and for |A • Uj\ 2 . 
Hence, the only relevant expressions are of the form ze 



-ik jy + 



re i(ip-k j y + ) an( j 



ze 



ikjy+ 



re 



-i(<p-kjy + ) 



When k\ = k 2 = k the energy density depends only on 
the combination (jp — ky + ) and r showing that the only effect of the dependence 
on time is a rotation and, in consequence, the independence of the energy per unit 
length on y + (or x° for given x 3 ). The other cases guaranteeing this are a\ — 
and a 3 = 0. 



4.2.3 The anti-holomorphic solution 

Finally we look at the corresponding anti-holomorphic solution. Such a solution 
derived from the system ( 4.34[ > takes the form 

Ul ^' y+) ~ p2f 3 ai a 3 + ze ik ^{2a 3 + a 2 ze ik ^+) 1 ' ' 

Plh a4 e^ + (a 1 + 2^e^+) 

U2[z,y+) - p2 j^ - aia3 + z e ik iy+ (2a 3 + a 2 ze ik *y+y { } 

Note that, like for the 'mixed case', we have to require that a 2 7^ as otherwise 

n U = 0, = ^M \ k } ala2 ' 

(ai + aj) 2 

In the next subsection we will produce an explicit example of these field con- 
figurations and discuss some of their properties. To avoid the problems mentioned 



13 



above our example will have a 2 ^ 0. Note that in such a case the conditions of the 
independence of the energy per unit length on y + are the same as for the mixed 
solution. 



4.2.4 An example 



In our example we start with the set of functions (4.34) for which we have chosen 
the following values of parameters: a\ = 2.5, a 2 = 0.6, 03 = 1.0, 04 = 0.01 
ki = 1.0 and k 2 = 2.0. The topological charge of the holomorphic solution is then 
Qtop = 2. The topological charge density at x° = (and for x 3 = 0) has two peaks 
- one of them is localized at z = 0, the other a bit further out - see Fig. [7} For the 
holomorphic solution the topological charge density is proportional to the energy 
density and this leads to the energy per unit length being given by SuM 2 !^. 
The integrand T-L^/SM 2 is sketched in Fig. [8] The contribution coming from 
the waves /8M 2 is plotted in Fig. |9| The mixed solution generated by the 



application of the P z operator according to (3.8) leads to a solution which has 
Qt op = 2 - 2 = and X (1) =2 + 2 = 4. As is easy to see from Fig. [7] the 
application of P z has changed two holomorphic peaks into two anti-peaks and in 
addition it has generated two new peaks. The energy density H^ 1 ' thus has four 
peaks and W- 2 ' only three (with the zero in the place of the fourth one). 

The next application of the P z operator changes two peaks of the topological 
charge density into two anti-peaks and annihilates the previous anti-peaks. Thus 
the anti-holomorphic solution is characterized by Q top = —2 and ZW = 2. The 
contribution to the energy per unit length 8M 2 X^ is the same as for the initial 
(holomorphic) case. Nevertheless, the total energies per unit length for these two 

solutions differ since for solutions of the CP 2 model the integrals are different 

(2) (2) 

(Z^ Q i 7^ ^anti-hoi)- Let us note that our case has a time-dependent energy per 
unit length (calculated by the integration over the x l x 2 plane). It implies that the 
dependence of the energy density on y + is highly nontrivial. However, the energy 
per unit length is a periodic function of y + . Only for some special cases, like 
k\ = k 2 etc. the energy per unit length is constant and so does not depend on y + . 
The time dependence of the energy density for all three solutions is shown in Fig. 



10 The energy density for the mixed solution for x 3 = and xq = it/ 4, xq = ir, 



xq = 77r/4 is plotted in Fig. 11 For the case x° = n the peaks are maximally 
separated (this is not very clear without a detailed study of some other values x°). 



In this case the energy takes its maximal value, see Fig. 10 
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Figure 1: The tube solution. The part (n 1 , n 2 ) (left) and the component n 3 (right) 
of the isovector n for a\ = 0, a 2 = 0, a 3 = 2, x° = 0, x 3 = and k 2 = 2. The 
minimal value n 3 = — 1 occurs at the point x 1 = and x 2 = 0. 




4 r 1 ^- 4 

Figure 2: The energy density of the tube solution - the topological part (left) and 
the wave part (right). Here a x = 0, a 2 = 0, a 3 = 2, a; = 0, a; 3 = and k 2 = 2. 
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Figure 3: The spiral solution. The part (n 1 , n 2 ) (left) and the component n 3 (right) 
of the isovector n for a% = 2, a 2 = 1, 0.3 = 0, x° = 0, x 3 = and k\ = 1. 
The minimal value n 3 = — 1 occurs at the point rr 1 = — a\ cos (kiy + ) and x 2 = 
— ai sin (kiy + ); here x 1 = —2, x 2 = 0. 




Figure 4: The energy density of the spiral solution - the topological part (left) and 
the wave part (right). Here a\ = 2, a 2 = 1, a 3 = 0, and x° = 0, x 3 = and 
ki — 1. The maxima of the energy density for both contributions are located at 
the same point on the plane x 1 x 2 corresponding with the minimum of n 3 (see Fig 
[3]); here x 1 = -2, x 2 = 0. 
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Figure 5: The CP 1 solution with all 7^ 0. The part [n 1 ^ 2 ) (left) and the 
component n 3 (right) of isovector n for a\ = 2, a 2 = 1, 03 = 3, x° = 37r/4, 
x 3 = 0, ki = 1 and A; 2 = 2. 




Figure 6: The energy density of the CP 1 solution with all ^ - the topological 
part (left) and the wave part (right). Here a\ = 2, a 2 = 1, a 3 = 3, x° = 37r/4, 
a; 3 = 0, k x = 1 and A; 2 = 2. 
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Figure 7: The functions 7rp top (where p top is a topological charge density) for x = 
0, x 3 = Q,ai = 2.5, a 2 = 0.6, a 3 = 1.0, a 4 = 0.01, fc a = 1 and k 2 = 2. The left 
picture corresponds to the holomorphic solution, the central picture corresponds 
to the mixed solution and the right picture to the anti-holomorphic solution. 




Figure 8: The functions /8M 2 (proportional to topological contribution to the 
energy density) for x° = 0, x 3 = 0, a\ = 2.5, a 2 = 0.6, a 3 = 1.0, a 4 = 0.01, 
ki — 1 and k 2 = 2. The left picture corresponds to the holomorphic solution, 
the central picture corresponds to the mixed solution and the right picture to the 
anti-holomorphic solution. 
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Figure 9: The functions "H^ 2 '/8M 2 (proportional to wave contribution to the en- 
ergy density) for x° = 0, x 3 = 0, ax = 2.5, a 2 = 0.6, a 3 = 1.0, a 4 = 0.01, 
k\ = 1 and &2 = 2. The left picture corresponds to the holomorphic solution, 
the central picture corresponds to the mixed solution and the right picture to the 
anti-holomorphic solution. 




Figure 10: The integral as the function of x° £ [0, 2n]. The other parameters 
read: x 3 — 0, a\ — 2.5, a 2 = 0.6, a 3 = 1.0, a 4 = 0.01, ki — 1 and A; 2 = 
2. The left picture corresponds to the holomorphic solution, the central picture 
corresponds to the mixed solution and the right picture to the anti-holomorphic 
solution. 
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Figure 11: Time evolution of the mixed solution for a\ = 2.5, a 2 = 0.6, a 3 = 
1.0, a 4 = 0.01, ki = 1 and k 2 = 2. The functions H (1) /8M 2 (left column) 
and 'US 2 ' /%M 2 (right column) have been considered for x 3 = at the moments 
x° = 7r/4 (first row), x° = n (second row), a; = 77r/4 (third row). Their values 
for x° = are sketched at the central pictures of Fig [8] and Fig [9} 
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5 Conclusions and Further Comments 



In this paper we have demonstrated that the CP model in (3+1) dimensions has 
many classical solutions. Our construction has been based on the observation that 
one can generalise ideas used in the construction of solutions of the CP N model in 
(2+0) dimensions and generate vortex and vortex-antivortex like solutions of this 
model in (3+1) dimensions. Like for the model in (2+0) dimensions we can gen- 
erate these solutions from field configurations described by polynomial functions 
of x 1 + ieix 2 . This time the coefficients of these functions could be also functions 
of x 3 + e 2 x . The energy of such configurations is infinite (as the energy density 
is independent of a; 3 ) and so we interpret these solutions as describing systems of 
vortices and antivortices. 

Of course our expressions solve equations in (3+1) dimensions and they also 
determine the dynamics of these vortices. 

In this paper we have only looked at the simplest solutions (corresponding 
to very few vortices) with the time dependence being described by simple phase 
factors. Even in this case the observed dynamics is quite complicated and has 
exhibited various interesting properties. In particular, we have shown that the 
vortices can rotate in space (physical and internal) and their energy per unit length 
of the vortex can vary in time. During this time evolution some vortices can 
shrink to delta functions and then expand again often being characterised by a 
very periodical behaviour. 

One other unusual property is their dependence on the distance between the 
vortices: the energy density of two vortices can depend on the distance between 
them and can possess a minimum at a specific value of this distance. This sug- 
gests that vortices which are located at non-minimal distances may be unstable 
and so could try to reduce their energy per unit length by moving towards this op- 
timal configurations. However, their configurations are solutions for any distance 
as their infinite 'inertia' stops them from moving towards each other without an 
external push. 

We are now looking at other properties of these and other solutions. 
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